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Abstract 

Purpose The results of published Life Cycle Assessments 
(LCAs) of biofuels are characterized by a large variability, 
arising from the diversity of both biofuel chains and the 
methodologies used to estimate inventory data. Here, we 
suggest that the best option to maximize the accuracy of 
biofuel LCA is to produce local results taking into account 
the local soil, climatic and agricultural management factors. 
Methods We focused on a case study involving the produc¬ 
tion of first-generation ethanol from sugar beet in the Pic¬ 
ardy region in Northern France. To account for local factors, 
we first defined three climatic patterns according to rainfall 
from a 20-year series of weather data. We subsequently 
defined two crop rotations with sugar beet as a break crop, 
corresponding to current practice and an optimized manage¬ 
ment scenario, respectively. The six combinations of climate 
types and rotations were run with the process-based model 


Responsible editor: Seungdo Kim 

C. Bessou (E) 

CIRAD, 

Performance of Perennial Cropping Systems Research Unit, 
Avenue Agropolis, 

34398 Montpellier, France 
e-mail: cecile.bessou@cirad.fr 

S. Lehuger 

INRA, Environment and Agricultural Crop Research Unit, 
78850 Thiverval-Grignon, France 

B. Gabrielle 
AgroParisTech, INRA, 

Environment and Agricultural Crop Research Unit, 

78850 Thiverval-Grignon, France 

B. Mary 

INRA, US 1158 Agro-Impact, 

02007 Laon-Mons, France 


CERES-EGC to estimate crop yields and environmental 
emissions. We completed the data inventory and compiled 
the impact assessments using Simapro v.7.1 and Ecoinvent 
database v2.0. 

Results Overall, sugar beet ethanol had lower impacts than 
gasoline for the abiotic depletion, global warming, ozone layer 
depletion and photochemical oxidation categories. In particu¬ 
lar, it emitted between 28 % and 42 % less greenhouse gases 
than gasoline. Conversely, sugar beet ethanol had higher 
impacts than gasoline for acidification and eutrophication 
due to losses of reactive nitrogen in the arable field. Thus, 
LCA results were highly sensitive to changes in local con¬ 
ditions and management factors. As a result, an average 
impact figures for a given biofuel chain at regional or national 
scales may only be indicative within a large uncertainty band. 
Conclusions Although the crop model made it possible to 
take local factors into account in the life-cycle inventory, best 
management practices that achieved high yields while reduc¬ 
ing environmental impacts could not be identified. Further 
modelling developments are necessary to better account for 
the effects of management practices, in particular regarding 
the benefits of fertiliser incorporation into the topsoil in terms 
of nitrogen losses abatement. Supplementary data and model¬ 
ling developments also are needed to better estimate the 
emissions of pesticides and heavy metals in the field. 

Keywords Agricultural practices - Biofuel • CERES-EGC • 
Ethanol • Greenhouse gases • Local LCA • N 2 0 • NOE2 • 
Process-based model • Sugar beet 

1 Introduction 

The development of renewable energy sources has recently 
been fostered and promoted as an alternative to the use of 
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fossil resources to mitigate greenhouse gas (GHG) emis¬ 
sions and climate change worldwide (IPCC 2011). While 
no single renewable energy source may claim to be the 
perfect solution, it is their best combination in relation to 
local resources that may lead to a carbon neutral energy mix. 
Among these renewable sources, biofuels have been largely 
promoted worldwide. In Europe, the renewable directive 
(2009/28/EC, European Commission 2009) aims at substi¬ 
tuting 10 % of its fossil fuel consumption for transport with 
biofuels by 2020. In order to be eligible for this share and to 
benefit from fiscal incentives, biofuels should comply with a 
set of sustainability criteria, including minimum GHG sav¬ 
ings of 35 % compared to equivalent fossil fuels. These 
criteria were imposed to ensure a net benefit in terms of 
environmental impacts with biofuels, following the contro¬ 
versy on their life-cycle GHG emissions (Quirin et al. 2004; 
Farrell et al. 2006). 

Life Cycle Assessment (LCA) is the most commonly 
used framework to assess the environmental impacts of 
products throughout their life cycle, since its holistic 
approach makes it possible to avoid pollution trade-offs across 
the life-cycle stages or impacts. However, the complexity of 
biofuel chains has evidenced some limits of this methodology, 
in particular because system boundaries are variable across 
studies, and because the assumption that emissions are pro¬ 
portional to inputs is often violated for field emissions. These 
emissions result from complex biogeochemical processes that 
are strongly dependent on local soil, climatic and management 
factors. For instance, the emissions of nitrous oxide (N 2 0), a 
gas with a global warming potential (GWP) 298 times larger 
than C0 2 (Forster et al. 2007), depend on a number of the 
above-mentioned factors and are therefore highly variable at 
the plot scale (Parkin 1987; Mosier et al. 1996; Farquharson 
and Baldock 2008). As a consequence, the linear regression 
model used in the IPCC Tier 1 (2006) estimation method, 
which only accounts for soil N input rates should only be 
expected to provide an order of magnitude for these emis¬ 
sions, associated with a wide uncertainty band (-70 % 
to +300 %; IPCC 2006). 

Also, many studies on biofuel chains are focused on spe¬ 
cific impacts (e.g., global warming) rather than the range of 
impacts that may be potentially addressed by LCAs (Quirin et 
al. 2004; von Blottnitz and Curran 2007), thus conveying an 
incomplete picture of biofuels sustainability. Finally, given the 
expanding demand for biofuels and the limited available land 
area, concerns have risen that biofuels may lead to land use 
changes, whose impacts are not yet fully assessed within the 
classical (attributional 1 ) LCA framework (Reinhard and Zah 

1 An attributional LCA is a descriptive LCA assessing the impacts of 
an added produced Functional Unit=FU (status-quo+FU) without 
considering the changes induced by this added FU as it is done in a 
consequential or change-oriented LCA (FU+consequential modifica¬ 
tions of the initial state due to the added FU). 


2009). Because of all these limitations of LCA, confusion and 
controversy arose on the environmental interest of biofuels 
(Crutzen et al. 2008; Fargione et al. 2008; Bessou et al. 2010a). 

Assessing the environmental impacts of a biofuel chain 
with enough accuracy so that a decision can be made with 
respect to the above-mentioned sustainability criteria requires 
(1) the definition of common system boundaries and method¬ 
ological choices and (2) a reduction of the uncertainty of 
inventory data, which may be achieved by more accurate 
emission models. There are currently ongoing international 
projects seeking to harmonize the LCA methodology 
(PAS2050 by BSI 2008; ISO 14067 2 ; ILCD by European 
Commission-JRC-IES 2011), while past studies have dealt 
with system boundaries for biofuel chains (Quirin et al. 
2004; Farrell et al. 2006). Other studies also used agro¬ 
ecosystem models to generate input data for the life-cycle 
inventory, thereby reducing the uncertainty of emissions from 
agricultural fields (Gabrielle and Gagnaire 2008; Langevin et 
al. 2010). Nevertheless, these model outputs were either spe¬ 
cific to a set of particular sites (for the latter studies) or related 
to country averages in terms of management data (JRC 2008; 
Smeets et al. 2009). When dealing with biofuels, data repre¬ 
sentative from the feedstock supply area should be sought. 

Our objective was hence to produce such a 'local LCA' 
for first generation ethanol from sugar beet, and to test the 
sensitivity of its results to changes in local environmental or 
management factors. We used the Picardy region in Northern 
France as a case study, involving sugar beet as a feedstock. 
Sugar beet has been produced on a large scale for decades in 
this area. To account for in the effect of local factors, we 
simulated crop growth and N emissions to the environment 
with the dynamic agro-ecosystem model CERES-EGC 
(Gabrielle et al. 2006). This works addresses regional varia¬ 
tions and data quality but does not address inconsistency in 
system boundary and the linearity assumptions imposed on 
the parameters not modelled. 


2 Material and methods 

2.1 Goal and system boundary 

The aim of this study was to conduct a full-blown attribu¬ 
tional LCA of the production and combustion of first- 
generation ethanol from sugar beet. We compared ethanol 
from sugar beet with its fossil-based equivalent, gasoline, to 
examine the benefits of substituting the latter with bio¬ 
ethanol, however strictly in an attributional framework. 
System boundaries were Well-to-Tank, i.e., including all 
flows of resources, energy and pollutants from the extrac¬ 
tion of raw materials up to ethanol dehydration. The chosen 

2 ISO 14067 still under development. 
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functional unit was 1 MJ of biofuel energy content, using 
the low heating value (LHV) of ethanol. This was calculated 
from the potential energy that would be delivered by the fuel 
at the refinery gate if its combustion in the car engine was 
complete. Since we focused here on the impacts of biofuel 
production and more specifically on its agricultural feed¬ 
stock, we did not include the steps of fuel blending and 
distribution to the filling station, as both would be similar 
between the two fuels. Energy embedded in farm machinery 
and capital equipment was included. On the contrary, the 
caloric intake and transportation of farm workers were not 
included (Farrell et al. 2006). Missing reliable inventory 
data prevented a complete inventory of some packaging, 
enzyme and yeast inputs for the fermentation processes. In 
terms of land use, we considered that sugar beet was grown 
on existing crop land, so that no land-use changes occurred 
(PAS2050 by BSI 2008). From an attributional perspective, 
we assessed the impacts of producing sugar beet ethanol 
without accounting for any possible displacement of other 
crops to compensate for the resulting reduction in sugar 
production. This approach was in accordance with the un¬ 
derlying assumption that there is an overproduction of sugar 
in this region, as a result of the sugar reform in the European 
Union. We did not account for biogenic C0 2 fixation and 
emissions since no emissions from land use change were 
considered. The system boundaries for the two fuel chains 
are presented (Fig. 1). 

The system was extended to take into account savings 
due to the substitution by co-products as recommended by 
LCA-ISO standards (ISO 14044:2006). Over the year, sugar 
beet ethanol is produced from different raw materials: sugar 
beet juice and fermentation slurry during the campaign (in 


the 3 months following the harvest of sugar beet crops), and 
syrup and slurry during the remaining time interval (inter¬ 
campaign). The inputs and outputs of the refinery-distillery 
processes were calculated based on the mass balances of the 
processes from raw material to the various final products 
over the year. The co-products of sugar beet ethanol con¬ 
sidered included pulp, lime-carbonation residue, vinasse and 
sugar. For each of these co-products, alternative chains were 
designed according to local practices and the expertise of 
professionals from the region. The avoided production of 
mineral fertilisers was calculated based on the chemical 
composition of co-products and their mineralization rates, 
depending on the site and time of spreading (Decloux et al. 
2002; ITB 2002). The equivalent of pulp was silage maize, 
given the close dry matter contents (around 25 % DM) and 
low protein values of the two products that leads farmers to 
use them indifferently in feeding rations (Comite National 
des Co-Produits/ADEME 2001). In state of the art of sugar 
beet sugar-and-ethanol production plants in France, sugar 
represents the main output. For the allocation between sugar 
and ethanol, we intrinsically used a mass-based procedure 
through the mass balance calculation. It is justified by the 
fact that the only locally available substitute for sugar 
involved the same sugar beet chain. Given this case, the 
use of system expansion and from mass-based allocation 
would give the same results. 

2.2 Temporal and geographical scopes 

LCA typically provides a ’snapshot’ of the studied systems at 
a given point in time. Here, the reference period for the data 
related to industrial processes was the 2000-2010 time 


Fig. 1 Well-to-tank system 
boundaries of compared fuel 
chains 
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interval, i.e., involving state-of-the-art technologies current¬ 
ly in use at commercial level. No prospective scenarios 
involving possible near-term improvements over these tech¬ 
nologies in terms of energy efficiency or pollution abate¬ 
ment were analysed. In particular, we did not consider any 
possibility of capturing fermentation-derived C0 2 . The ref¬ 
erence period for agricultural productions was the 1988- 
2009 time interval to account for sufficient climate variabil¬ 
ity in the modelling work. Crop management data were 
representative of a more recent period running from 2000 
to 2010. 

The agricultural production site considered was the 
INRA Agro-Impact experimental station at Estrees-Mons 
in the Picardy region (49°80'N, 3°60'E), where the mean 
annual precipitation and air temperature are 667 mm and 
9.6 °C, respectively (Guerif et al. 1997). The size of the 
experimental site is 170 ha. The annual total sugar beet area 
represents 15 ha. The soil is an Orthic Luvisol (FAO 1998), 
which represents 20 % of the arable soils in this region 
(Gagnaire et al. 2006). In order to better capture the emis¬ 
sion variability due to the local production factors, climatic 
and agricultural management data were collected from this 
station over 21 years. In this area, sugar beet is a common 
break crop rotated with wheat and maize. None of these 
crops are usually irrigated. 

2.3 Life cycle inventory of the agricultural production 
2.3.1 Sugar beet cropping system 

Field emissions are linked to physicochemical and microbial 
transformations in soils. They can be classified as direct or 
indirect according to their occurrence in the causal chain 
from the application of agricultural inputs to the final 
impacts incurred. For instance, primary emissions of 
leached nitrate are defined as direct field emissions, whereas 
N 2 0 emissions ensuing from this primary emitted nitrate but 
occurring beyond the agricultural field itself are considered 
indirect. Field emissions are related to reactive compounds 
either applied to the field as fertilisers and pesticides, or 
originating from the mineralization of crop residues. Most 
field-applied fertilisers and pesticides are managed over one 
crop cycle. However, some fertilisers such as lime may be 
managed at the crop rotation level. In our case study, lime 
was managed over a four-crop rotation, hence one fourth of 
the input rate and spreading costs were allocated to sugar 
beet. 

In order to account for all field emissions, we considered 
sugar beet as part of the sugar beet-winter wheat-maize- 
winter wheat rotation, which is common in the area. This 
rotation was implemented and monitored at the study site for 
more than 20 years. We simulated this rotation with 
CERES-EGC during the 1988-2009 period, meaning that 


each run consisted of five repetitions of the rotation. To 
obtain an occurrence of sugar beet for all individual years 
of the 1988-2009 period, we shifted the starting year from 
1988 to 1991, resulting in four model runs. For all crops, 
residues were incorporated by deep ploughing. Emissions 
due to residue mineralization and nitrate leaching were 
allocated to each crop of the rotation by considering that N 
mineralization was the main source of inorganic N in 
between two growing seasons. We first defined the time 
interval for this allocation as the period running from the 
sowing of the crop of interest to the first fertiliser application 
on the following crop. In terms of tillage, each crop was thus 
allocated the incorporation of its harvest residues. Since the 
studied rotation included some 8 months of bare soil be¬ 
tween wheat harvest and the sowing of sugar beet or maize, 
we allocated part of this bare soil period to these spring 
crops. This discrimination resulted in a final splitting of the 
whole rotation into four 1-year periods, starting and stop¬ 
ping on February 15 each year. N inputs were calculated 
from a balance of crop requirements and supply from soil, 
including residue returns from the preceding crops. Hence, 
the green residues from sugar beets left in the field after 
harvest were not considered as a co-product of ethanol in the 
LCA. 

The current crop rotation is not optimal in terms of soil 
cover. Introducing a cover crop between winter wheat and 
sugar beet, for instance, may reduce the bare soil period by 
some 85 days (ITB 2002). This cover crop is thereby 
expected to reduce N losses by taking up N prior to the 
winter drainage period, and to make it available to the 
following crop after the cover crop has been ploughed in. 
The combined effect of cover crop on reducing nitrate 
leaching and increasing efficiency of N uptake by the crops 
can be significant over the years (Constantin et al. 2010). N 
losses may be further reduced by adjusting fertiliser N 
management. In particular, the soil incorporation of the first 
fertiliser N application together with the seeds can save as 
much as 15 kg N ha -1 year -1 (ITB 2002). Beside our 
baseline scenario of the current rotation, we hence tested 
an 'optimized rotation', consisting of two aspects: (1) the 
introduction of a cover crop within the rotation between 
winter wheat and sugar beet or maize and (2) the incorpo¬ 
ration of N inputs upon sowing allowing the above men¬ 
tioned savings of 15 kg N ha -1 year -1 in both sugar beet and 
maize fields. In the study, we extended the cover crop, i.e., 
mustard, from August 15 until January 5, i.e., 144 days, so 
that it would be killed by the frost without requiring herbi¬ 
cides. The consequences for sugar beet were accounted for 
through the resulting decrease of N losses between February 
15th and the sowing in April as simulated by the crop 
model. The costs of cover crop establishment, i.e., sowing 
and mechanical destruction were equally split between sugar 
beet, winter wheat and maize. Annual P and K fertilization 
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rates in sugar beet plots were 70 kg ha -1 P 2 O 5 , and 
110 kg ha -1 K 2 0. These inputs, as well as pesticides and 
tillage operations (except for the introduction of the cover 
crop) did not vary between the two management scenar¬ 
ios. N fertiliser rates and modelled yields are shown in 
Table 1. 

2.3.2 Modelling of field emissions 

To take into account local soil and climatic factors in the 
LCA, we modelled feedstock production with the agro¬ 
ecosystem model CERES-EGC (Gabrielle et al. 2006). This 
process-based model simulates crop growth by considering 
the type of soil, climate conditions, plant variety, and man¬ 
agement practices at the field scale. It also predicts sub¬ 
stance losses outside the ecosystem such as nitrate leaching, 
ammonium volatilization (Genermont and Cellier 1997), 
and nitrous oxide emissions through nitrification and deni¬ 
trification which are predicted with the NOE sub-model 
(Henault et al. 2005). In order to improve its accuracy, we 
implemented the NOE2 sub-model (Bessou et al. 2010b) in 
CERES-EGC and proceeded to a Bayesian calibration of the 
22 parameters of NOE2. This calibration makes it possible 
to reduce significantly the spread in parameter distribution 
and the subsequent uncertainty on N 2 0 emissions by updat¬ 
ing these distributions against a priori probability distribu¬ 
tions of parameter values gathered from other sites (Lehuger 
et al. 2009). This calibration reduced the root-mean square 
Error on daily and annual fluxes by 15 % and 63 %, respec¬ 
tively, when compared to the uncalibrated simulations 
(Fig. 2), based on experimental data collected in our study 
site over sugar beet (Bessou et al. 2010b). Calibrated param¬ 
eters were then used for the simulations of sugar beet crop¬ 
ping system with CERES-EGC (Leviel 2000; Leviel et al. 
2003). Despite our efforts, CERES-EGC still underesti¬ 
mated emissions in 2008-2009 due to several combined 
factors (Bessou et al. 2010b). Other sources also showed a 
marginal error of roughly 20 % with CERES-EGC (Henault 
et al. 2005; Gabrielle et al. 2011). This error remains lower 
than the uncertainty associated with IPCC factors. 

We used daily climate data from a weather station located 
on site, and input soil and management data for the 21 -year 
simulation period (1988-2009). According to climate 
records over this period and annual rainfall, we identified 
three climate patterns: (1) median years (R50; there were 
9 years in this class in the simulation period); (2) inferior 
quartile years, ’drier’ than median years (R25; 5 years); and 
(3) superior quartile years, ’wetter’ than median years (R75; 
5 years) (Fig. 3). Average values (arithmetic means) of 
modelled yields and emissions were calculated for the six 
combinations of rotation management and climatic years, 
noted R25, R50, and R75 for the baseline crop rotation, and 
R25 optim, R50 optim, and R75 optim for the optimized 


crop rotation. The resulting data are given in Table 1, and 
were fed into the life cycle inventories. An overview of 
the study approach and tested parameters is given in 
Fig. 4. 

2.3.3 Field emissions not modelled with CERES-EGC 

Changes in the spatial and temporal scales between direct 
and indirect field emissions make it particularly difficult to 
assess indirect emissions. Mechanisms are poorly character¬ 
ized and indirect emissions are therefore not modelled by 
CERES-EGC. We thus used IPCC coefficients as an alter¬ 
native to estimated indirect emissions, and also to estimate 
C0 2 emissions due to lime application. 

Methane emissions are not considered by CERES-EGC, 
but were ignored since the soil was well drained and unlike¬ 
ly to be a net emitter or consumer of CH 4 (BIOIS 2008). 
Finally, we accounted for heavy metals brought to the soil 
by fertilisers following the corresponding inventory method 
of Ecoinvent. Due to the lack of reliable data on the heavy 
metals brought to the soil in pesticides or exported with the 
harvested biomass, we decided not to account for them. The 
modelling of direct and indirect emissions linked to 
agricultural productions in this study in summarized in 
Table 2, also showing the differences with their respec- 
five models in the Ecoinvent documentation (Nemecek 
and Kagi 2007). 

2.4 Life cycle inventory of the industrial phase 

Background data on primary energy sources, infrastructures, 
machines, industrial production and disposal were mainly 
taken from Ecoinvent v2.0 databases (Swiss Centre for 
LCI). We used the data provided on the average French mix 
for electricity from the grid, as well as European average data 
for crude oil supply, oil refinery and energy embedded in 
infrastructures and machines. Transports of inputs to the agri¬ 
cultural field encompassed that of co-products from the refin¬ 
ery to the field, as well as saved transports for the products 
they substituted, except for maize silage which was produced 
on farm. Transport distances corresponded to current practices 
at the study site: 65 km for the supply of agricultural inputs to 
the farm, 22 km for the transport of sugar beet roots to the 
sugar refinery, and 200 km for the other inputs of the refinery. 
Conversion data was taken from the Arcis-sur-Aube sugar 
refinery-distillery, which is located in the Champagne- 
Ardennes region neighbouring Picardy. This refinery-distillery 
is representative of the state of the art of such production units 
in France (ADEME 2010). 

3 The Ecoinvent Centre—a Swiss Competence Centre of ETHZ, 
EPFL, PSI, Empa and ART—is the world's leading supplier of consis¬ 
tent and transparent life cycle inventory (LCI) data of known quality 
with the database ecoinvent data v2.2. See http://www.ecoinvent.ch/. 
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2.5 Impact characterization 

We used the CML 2000 Baseline v2 method (Guinee et al. 
2001) for the impact characterization, similarly to the refer¬ 
ence LCA study commissioned by the French Environment 
Agency (ADEME 2010). CML 2000 encompasses the 
impact categories listed in Table 3. Impact indicators were 
not subsequently normalized or weighted since these steps 
are partly subjective and uncertain, and hamper the interpre¬ 
tation of the LCA results (Guinee et al. 2001). Weighting in 
particular is not allowed by LCA-ISO standards when pub¬ 
licly comparing products, which is our objective here. LCA 
calculations were carried out with the SimaPro 7.1 software 
(PRe Consultants bv, NL). 


3 Results and discussion 

3.1 Biofuel versus fossil fuel 

Figure 5 compares the environmental impacts of the pro¬ 
duction and the pseudo complete combustion of sugar beet 
ethanol and gasoline (low sulphur grade) on a MJ basis. Six 
scenarios of sugar beet ethanol productions were compared, 
corresponding to the six combinations of climate patterns 
(3) and crop rotations managements (2). Overall, the 
impacts of sugar beet ethanol were lower than those of 
gasoline in four categories: abiotic depletion, global warm¬ 
ing, ozone layer depletion, and photochemical oxidation. 
Conversely, gasoline outranked bio-ethanol for acidifica¬ 
tion, eutrophication, human toxicity, and fresh water or 
terrestrial ecotoxicity. Feedstock production contributed 
more than 75 % of total impact indicators for acidification 
and eutrophication for ethanol, with emissions of reactive N 
from fertiliser inputs in the field being the top contributors. 
Despite substantial reduction in nitrate leaching with the 
optimized rotation (see Table 1), field emissions dominated 
the eutrophication impact. Reactive N also contributed to a 
large extent to GHG emissions in the field, while agricul¬ 
tural production contributed 25-68 % of the other impact 
categories. The industrial conversion processes contributed 
to a large extent (more than 55 %) to abiotic depletion, 
ozone layer depletion, aquatic ecotoxicity and photochemi¬ 
cal oxidation. Across impact categories, transport only con¬ 
tributed to 5-10 % of total impact indicators. 

Sugar beet ethanol produced in Picardy lead to savings of 
28^-2 % GHGs and 20-28 % fossil resources compared to 
gasoline, on a MJ basis. The mean GHG emissions across 
scenarios amounted to 5.78x10 kg C0 2 eq MJ , 
corresponding to a 35.7 % of GHG saving compared to 
gasoline. These emissions are larger than the France-wide 
average than those reported for sugar beet ethanol by 
ADEME, at 3.04x 10“ 2 kg C0 2 eq MJ -1 (ADEME 2010), 
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Fig. 2 Results of the Bayesian calibration of CERES-NOE2 (best set 
of parameter values), simulated fluxes in black, observed in grey a 
sugar beet plot in 2007: RMSE (N 2 0 daily emissions): 6.97 gN-N 2 0 
ha 1 day _1 ; RMSE (cumulated N 2 0 emissions): 138.89 gN-N 2 0 ha -1 ; 

but close to the European average value of 5.75x10 kg 
C0 2 eq MJ 1 published by JRC (2008), where sugar beet 
pulp was also used as cattle feed. Across our six feedstock 
scenarios, the global warming, eutrophication, and human 
toxicity indicators laid within -51 % and +113 % of the 
ADEME study (2010). Global warming impacts were all 
higher than ADEME (2010) (+80 % to +113 %), as well as 
human toxicity impacts (+24 % to +52 %), while eutrophi¬ 
cation impacts were either lower (-51 %) or higher (+23 %) 
depending on the feedstock scenario. However, both studies 
are not directly comparable due to differences in system 
boundaries and field emission models, and scope. The 
ADEME (2010) and JRC (2008) estimates apply to national 
to continental scales and are relevant to public policies at 
macro-economic level, while our results were more local 
and targeted at stakeholders such as farmers or chain oper¬ 
ators in the region. The sensitivity of LCA results to local 
factors highlights the relevance of this local LCA to define 
strategies at this level. 

In our study, the conversion of 1-ha sugar beet provided 
between 2.60 and 3.97 t ethanol across the scenarios, since 
ethanol only made up 35 % of the total outputs from the 
sugar refinery-distillery. Should this ratio rise to 100 %, as 
assumed in the JRC study (2008), ethanol yields per ha 
would be nearly twice as high (at 5.78 tha -1 ; ADEME/ 
DIREM 2002). Because the data collected from the sugar 
refinery-distillery were aggregated across the ethanol and 
sugar end-products and both chains were tightly entangled, 
it was impossible to test such a scenario and subsequently 
apply an energy-based allocation between ethanol and its 




b sugar beet plot in 2008-2009: RMSEP (N 2 0 daily emissions): 8.04 g 
N-N 2 0 ha 1 day -1 ; RMSEP (cumulated N 2 0 emissions): 266.78 gN- 
N 2 0 ha -1 


remaining co-products, as recommended in the European 
Renewable Directive (RED) methodology (2009/28/EC, 
European Commission 2009). At any rate, this scenario 
would not be representative of the current situation, since 
industrial processes have been up to now optimized to 
reduce energy costs for the combined production of sugar 
and ethanol during both campaign and inter-campaign peri¬ 
ods. Moreover, sugar has a higher added value than ethanol, 
making the 100 % ethanol unrealistic. To our understanding, 
this duality was not considered in the RED methodology. 
Despite its particular configuration, sugar beet ethanol in 
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Fig. 3 Rainfall patterns 
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industrial input 



Fig. 4 Overview of the study approach and tested parameters 


Table 2 Modelling of field natural emissions in our case study compared to modelling field emissions in Ecoinvent 


Emissions 
(kg ha 1 year -1 ) 

Modelled in our case study 

Modelled in Ecoinvent (Nemecek and Kagi 2007) 

Direct 

Modelled with CERES-EGC model 

Modelled with several operational models 

n-n 2 o 

Produced via denitrification and nitrification of ammonium 
nitrate as influenced by soil temperature, soil water content, 
soil compaction, substrate availability, and microbial potentials. 

1.25 % of N input (less volatilized N) in kg ha 1 


Source: Henault et al. 2005, Bessou et al. 2010 

Source: Schmid et al. 2000 


Modelled range: [0.8-1 % of N input] 

Single value 

N-NO x 

0.5 % of nitrified N 

N-NO x =21 % of N-N 2 0 


Source: Rolland et al. 2010 

Source: Schmid et al. 2000 


Modelled range: [0.5-0.6 % of N input] 

Single value 

N-NH 3 

Produced via volatilization as influenced by soil energy balance 

2 % of N input in kg ha 1 


Source: Gabrielle et al. 2006 

Source: Asman 1992 


Modelled range: [2.6-5.2 % of N input] 

Single value 

N-N0 3 - 

Through leaching as influenced by soil type and depth, soil 
water content, plant uptake 

24 % of N input in kg ha 1 for sugar beet 


Source: Gabrielle et al. 2006 

Source: Richner et al. 2006 


Modelled range: [3.7-31.2 % of N input] 

Single value 

Direct 

Not modelled with CERES-EGC model but with several 
operational models 

Modelled with several operational models 

C-C0 2 

From lime application 12 % of C-CaC0 3 year -1 with a large 
uncertainty (50 %), 

12 % being the maximum possible fraction 

Source: IPCC 2006 

Not available 

C-CH 4 

Not considered 

Not available 

P-P0 4 3- 

In kg P ha 1 year -1 

-through leaching in arable land (to groundwater): 0.07 

-through run-off in open arable land (to surface water): 0.175x 
(1+0.2/80) kg P 2 0 5 

-through erosion to rivers (not considered in our case study 
due to some lack of data) 

Source: Prasuhn 2006 


Indirect 

Not modelled with CERES-EGC model but with several 
operational models 

Modelled with several operational models 

N-N 2 0 

0.75 % of leached N0 3 “N+1 % of volatilized NO x -N and NH 3 -N 

2.5 % of leached N0 3 - -N+1 % of volatilized NH 3 -N 


Source: IPCC 2006 

Source: Schmid et al. 2000 
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Table 3 Baseline impact categories in CML 2000 v2.0 (Guinee et al. 2001) 

Impact categories 

Reference substances 

Main contributing substances 

Depletion of abiotic resources 

kgeq Sb 

Iron ore, Crude oil 

Acidification 

kgeq S0 2 

S0 2 , NO x , NH 3 

Eutrophication 

kgeq P0 4 3 “ 

P0 4 , N0 3 , most N compounds 

Global warming (GWP20, GWP 100) 

kgeq C0 2 

C0 2 , N 2 0 

Stratospheric ozone layer depletion 

kg eq CFC-11 

CFC 

Human toxicity 

kgeq 1,4-DB 

N 2 0, NO*, NH 3 , particles, S0 2 , acrylonitrile, cadmium 

Ecotoxicity: freshwater, marine, and terrestrial 

kgeq 1,4-DB 

Acrylonitrile, cadmium, particles, N 2 0, NH 3 , S0 2 

Photo-oxidant formation 

kgeq C 2 H 4 

CO, toluene, NO* 


Picardy would comply with the RED requirement of 35 % 
GHG savings. Further conversion of the current sugar co¬ 
product to ethanol would indeed lead to higher ethanol 
yields per hectare with identical field emissions, and conse¬ 
quently lower GHG emissions per MJ of ethanol. Neverthe¬ 
less, our results pointed out though that considering local 
factors, when data is available, might be crucial and more 
relevant than default parameters from the RED guidelines. 

Since assessing the life-cycle GHG emissions of biofuels 
is essential from a policy standpoint, we examined its sen¬ 
sitivity to the GWP of GHGs. We found an inconsistency 
between the value used in the JRC (2008) report and that 
reported in the RED e methodology description (Annex 
V.C), both cited by the same Directive. This was also 
reported in the Biograce calculator (IE 2011). These sources 
used GWP values (for a 100-year horizon) of 23 kg C0 2 eq 
per kg CH 4 and 296 kg C0 2 eq per kg N 2 0 (IPCC 2001) or 
25 kg C0 2 eq per kg CH 4 and 298 kg C0 2 eq per kg N 2 0 


(IPCC 2007). In our calculations, there were actually no 
discernible differences on the final GHG balances between 
the two sets of GWP values. Using GWPs for a 20-year 
period (considering higher concerns for immediate conse¬ 
quences and action needs) had a more visible effect. For the 
R50 scenarios, for instance, the largest differences occurred 
between the 20-year (2007) and 100-year (2001) perspec¬ 
tives, with GHG emissions of 6.20 x 10 -2 and 5.70 x 1(T 2 kg 
C0 2 eq per MJ, respectively. This difference would corre¬ 
spond to either 32.6 % or 36.8 % GHG savings compared to 
gasoline, respectively. We completed this sensitivity analy¬ 
sis by halving the C0 2 emission factor of applied lime, as 
suggested in the IPCC Guidelines (2006), but it only had a 
marginal effect on GHG budgets. Sugar beet development is 
sensitive to pH and to a risk of surface crust in the soil of 
Estrees-Mons, so that liming must be carefully managed. In 
our study site, the liming rates were minimal because of 
regular lime applications in the past. However, in the cases 


Fig. 5 Compared LCA 
results of sugar beet ethanol 
(6 combinations) and 
gasoline on a MJ basis 
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of soils with higher risks of surface crusting or insufficient 
liming record, rates could easily double and the consequen¬ 
ces on global warming impact become more important. This 
factor should be accounted for while identifying the best 
locations and managements for sugar beet farming. 

3.2 Impact of local production factors 

LCA results were sensitive to changes in climatic conditions 
during biomass production. As detailed in Table 1, local 
production factors influenced fresh matter yields, but also 
all N emissions that contribute to several impact categories. 
Consequently, local factors also had an influence on the 
final impacts. Compared to the sugar yields recorded on site 
over the 1990-2008 period, modelled yields were higher but 
with a larger standard deviation across years, but within the 
same order of magnitude as the observed yields. The ob¬ 
served yields as averaged over the drier, median and wetter 
years amounted to 78, 78 and 81 tha -1 (standard deviations: 
±11-13 %), compared to modelled yields of 76, 95, and 96 t 
ha -1 , respectively (standard deviations: ±15-44 %) (see 
Table 1). Climatic patterns affected impacts across catego¬ 
ries and management scenarios with a maximum relative 
difference of about 60 % in terms of GHG savings and an 
average difference per impact around 15 % (see Fig. 5). The 
largest difference occurred with eutrophication between the 
R25 and R50 climates. 

Drier years (R25) were more radically distinct from the 
other climates, and had the largest standard deviations, 
indicating a higher uncertainty in dry years. On a hectare 
basis (Fig. 6), the impacts of sugar beet ethanol under the 


R25 climatic pattern were lower than for the other patterns, 
whatever the impact category. Drier years resulted in lower 
direct emissions of all N compounds and subsequent indi¬ 
rect emissions, leading to reductions in the eutrophication, 
acidification, and global warming impacts. However, drier 
years (R25) were also characterized by lower biomass yields 
and ethanol production per hectare. Ethanol yields (kg) were 
20 % lower on average for the R25 scenario compared to 
R50 and R75, the latter showing similar ethanol yields per 
hectare (±2 %). On a MJ basis (see Fig. 5), R25 ethanol had 
therefore higher environmental impacts, except for eutrophi¬ 
cation, which was the most sensitive to rainfall patterns. In 
the drier years (R25), nitrate leaching was on average eight 
times lower than in the other years. Between the R50 and 
R75 scenarios, differences also appeared for the acidifica¬ 
tion and eutrophication impacts, which emphasize their 
sensitivity to climatic conditions. 

When comparing baseline and optimized cropping sys¬ 
tems on a MJ basis, the impact reductions between ethanol 
and gasoline varied from 10 % to 45 % across impacts (see 
Fig. 5). The largest difference occurred for eutrophication 
between R50 and R50 optim. The introduction of a cover 
crop reduced nitrate leaching more efficiently during medi¬ 
an years (R50 optim) than during drier (R25 optim) or 
wetter years (R75 optim). Across impact categories and 
scenarios, the optimized rotations had higher impacts, 
except for acidification (in the drier years) and eutrophica¬ 
tion (in the median to wet years). However, robust conclu¬ 
sions on the comparison of both rotations could not be 
drawn because of the trade-off between yield levels and 
field emissions, resulting in inconsistent rankings across 


Fig. 6 Compared LCA 
results of sugar beet ethanol 
(6 combinations) on a 
hectare basis 
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impact categories. On average, yields were 16 % lower with 
the optimized scenarios because of the reduced fertiliser 
inputs. The lower yields counter-balanced the lower field 
emissions occurring with the optimized system, which is a 
result of the strategic choices behind the optimization of the 
crop rotation. The latter implicitly considered land area as 
the main constraint for resource preservation: reducing fer¬ 
tiliser inputs per hectare aimed at abating N losses, while 
shortening the period under bare soil was expected to reduce 
nutrient losses. Yield optimization was therefore not a major 
factor in the management design although they certainly 
were in the LCA. This emphasizes the need to factor in 
ecological intensification targets when designing agricultur¬ 
al strategies, to avoid potential extensification and conse¬ 
quential increases in agricultural land area where land area is 
a limiting factor. This study, however, does not consider 
consequential impacts. 

The lower yields with optimized systems may be also 
partly attributed to a modelling artefact. Indeed, the relative 
benefit from topsoil fertiliser incorporation at sowing could 
not be modelled since CERES-EGC considers that fertiliser 
N is always incorporated into the topsoil. It is thus likely 
that the baseline scenario without incorporation should have 
higher nitrogen losses, in particular through ammonia vola¬ 
tilization. On the contrary, the impact of introducing a cover 
crop was better modelled, since nitrate leaching was reduced 
by more than 50 % in R50 and R75 compared to initial crop 
rotation (see Table 1). Differences in LCA results between 
the baseline and optimized rotations were hence the largest 
for the eutrophication indicator (see Figs. 5 and 6). Eutro¬ 
phication and acidification were the only categories for 
which optimized rotations performed better than baseline 
ones on both hectare and MJ bases, i.e., the reduction of 
impacts compensated for reduced yields. Across the three 
optimized scenarios, R50 optim had the least global impacts 
when considering both hectare and MJ bases. Despite the 
above shortcomings in the modelling of management effects 
of the optimized rotation, standard deviations were overall 
smaller for optimized scenarios than for the baseline ones. 
Moreover, score differences per category were larger for the 
baseline rotation than for the optimized one. Hence, the 
optimized rotation might help reduce risk of high environ¬ 
mental impacts in case of variable climatic conditions. 

Modelled field emissions were within the same order of 
magnitude as those calculated with the coefficients from 
IPCC Tier 1 or Ecoinvent, although generally slightly lower 
for N 2 0 and N0 3 _ , and larger for NH 3 and NO x emissions 
(see Table 1). Compared to our modelled results, IPCC Tier 
1 tended to overestimate these emissions, whereas Ecoin¬ 
vent modelling would lead to much lower figures. Using 
inventory data output by the CERES-EGC model made it 
possible to better assess the combined effects of local soil, 
climate and agricultural management factors than using 


default values per rotation scenario calculated with simpler 
methodologies. 

4 Conclusions 

The local LCA reported here provided detailed assessments of 
potential environmental impacts of ethanol from sugar beet at 
Estrees-Mons in the Picardy region in Northern France. 
Effects of local soil types, climatic conditions and manage¬ 
ment practices were visible for all impact categories and on 
both hectare and MJ bases. They emphasized the need to take 
into consideration these local factors when evaluating biofuel 
chains to provide guidance in their implementation. The 
environmental performance of biofuels was highly variable 
according to management and climatic conditions, within a 
span that could include threshold values such as the 35 % 
GHG binding savings of the RED directive for instance. The 
magnitude of these effects was particularly important for the 
eutrophication and acidification impacts, which are strongly 
dependent on field emissions and hence local conditions. 

GHG savings of sugar beet ethanol compared to gasoline 
ranged between 28 % and 42 % with a mean global warming 
indicator across scenarios of 0.0578 kg C0 2 eq per MJ, i.e., 
35.7 % of GHG savings. These 35.7 % savings might be 
critical because they are closed to the limit of 35 %, prereq¬ 
uisite of the European Directive (2009/28/EC), and lower 
than the 50 % limit to be enforced in 2017. In the EU 
Directive, biofuel chain performances were calculated with 
energy allocation for co-product handling, which might 
allow for more savings in the sugar beet case. Moreover, 
this indicator could vary largely according to the time-scale 
and GWPs chosen. 

The results on Ozone layer depletion left little scope for 
discussion here, since the method we used to characterize 
this impact (see Table 3) does not account for any contribu¬ 
tion of N 2 0 to stratospheric ozone depletion, while it is 
recognized that N 2 0 plays an important part in ozone chem¬ 
istry (Forster et al. 2007). This should be further investigat¬ 
ed in order to reflect better the impacts of the biofuel chains 
we tested here. We also remained cautious with the ecotox- 
icity impacts, since we could not account for the complete 
balance of heavy metals in our agro-ecosystems, or for a 
likely distribution of pesticides residues in the various envi¬ 
ronmental compartments. It stresses that despite the holistic 
nature of LCA, many environmental mechanisms are still 
not fully understood, nor clearly assessed with the current 
methodology (Mamy et al. 2010). 

Analysis on both hectare and MJ bases showed important 
trade-offs between low impact scenarios and low ethanol 
yields, which made it more complex to define optimal 
choices in terms of agricultural management. The ranking 
of the baseline and optimized crop rotations was inconsistent 
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across impacts as a result. It emphasizes the need to factor in 
ecological intensification targets when designing agricultural 
strategies, to avoid potential extensification and consequential 
increases in agricultural land area, especially in situations where 
land area is a limiting factor. LCAs may play a prominent role 
in designing such strategies. Climatic years with extreme rain¬ 
fall patterns (R25 and R75 scenarios) presented the worst 
impacts, while drier years led to more severe impacts on a MJ 
basis due to poor yields per hectare. R75 tended to increase 
several impacts such as eutrophication due to increased nitrate 
leaching, or acidification and global warming. This effect might 
be due to nitrogen cycling processes and microbial activities 
that are generally enhanced in wet conditions. 

Overall, the optimized scenarios had lower environmen¬ 
tal impacts than the other scenarios, although lower yields 
counteracted this positive effect to some extent. For eutro¬ 
phication in particular, optimized rotations would be very 
beneficial. If the risk of drought or heavy rainfalls were 
higher, optimized rotations would decrease the risks of high 
impacts, the indicators of these scenarios varying within a 
narrower range than those of baseline scenarios. 

Examining the sensitivity of LCA results to local factors 
makes it possible to define the probability distribution of 
environmental impacts around average results. This approach 
should provide guidance to decision-makers by identifying 
climatic and soil configurations that may lead to high envi¬ 
ronmental impacts. Further model developments would be 
also needed to improve the process-based modelling of the 
local factors and notably agricultural practices. In the current 
version of CERES-EGC, nitrogen fertilisers are distributed 
over the first 30 cm and nitrification and denitrification calcu¬ 
lated over the first 20 cm with average temperature and water 
content for two sub-layers across these 20 cm. The discretiza¬ 
tion of ammonium distribution and soil water and temperature 
within fine soil layers could improve the accuracy of model 
estimation (Bessou et al. 2010b). The modelling of N dynam¬ 
ics in soil also deserves further improvements by differentiat¬ 
ing between fertiliser incorporation depths and potentially 
improved sub-routine for N release from cover crop residues. 

In Picardy region, a new potential local factor is land use 
change. Indeed, experiments have been ongoing since 2006 
on comparing the agricultural potential of diverse biofuel 
feedstock from perennial crops (REGIX project). In the near 
future, some of these perennials, such as Miscanthus, may 
replace or displace sugar beet fields. 
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